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We give an overview of the scaling of density of quasi-particles and excess energy (heat) for nearly 
adiabatic dynamics near quantum critical points (QCPs). In particular we discuss both sudden 
quenches of small amplitude and slow sweeps across the QCP. We show close connection between 
universal scaling of these quantities with the scaling behavior of the fidelity susceptibility and its 
' generalizations. In particular we argue that the Kibble-Zurek scaling can be easily understood using 

, this concept. We discuss how these scalings can be derived within the adiabatic perturbation theory 

and how using this approach slow and fast quenches can be treated within the same framework. We 
' also describe modifications of these scalings for finite temperature quenches and emphasize the im- 

portant role of statistics of low-energy excitations. In the end we mention some connections between 
adiabatic dynamics near critical points with dynamics associated with space-time singularities in 
' the metrics, which naturally emerges in such areas as cosmology and string theory. 
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Continuous quantum phase transitions (QPTs) have been a subject of intense theoretical research in recent decades 
1-^3 ' (see e.g. Refs. [H-Q for overview). Unlike usual phase transitions driven by temperature, QPTs are driven entirely 
by quantum fluctuations. They are believed to occur in many situations as described later in this book. Quite 
recently a second order QPT was observed in a cold atom system of interacting bosons in an optical lattice. There 
a system of interacting bosons was driven in real time from the supcrfluid to the insulating phase 4], confirming an 
earlier theoretical prediction 5]. Up till now, Ref. 4] has provided probably the cleanest experimental confirmation 
of a QPT. The unifying property of all continuous (second order) phase transitions is the emergent universality and 
^ ■ scale invariance of the long-distance low energy properties of the system near the quantum critical point (QCP) [IJ]. 
This universality implies that low-energy properties of the system can be described by very few parameters, like the 
correlation length or the gap, which typically have power-law scaling with the tuning parameter characterized by 
critical exponents. These exponents are not sensitive to the microscopic details of the Hamiltonian describing the 
, system, but rather depend only on the so-called universality class to which a given phase transition belongs 
(~| • Recent experimental progress in preparing and manipulating out-of-equilibrium nearly isolated systems has stim- 
O [ ulated intense theoretical research on quantum dynamics in closed systems. In particular, such issues as sudden 
i quantum quenches in low dimensional systems adiabatic dynamics near QCPs [1, and connection of dynam- 

ics and thermodynamics in quantum systems jlOl . came to the forefront of theoretical research. A very natural 
^f) . question can be posed about non-equilibrium behavior of systems near QPTs. Since in equilibrium second order 
^ ' phase transitions are characterized by universality, one can expect also universal behavior in driven systems in the 
. vicinity of a QCP. Such universality can be expected, for example, if the system near the QCP is a subject to small 
amplitude, low frequency modulation of an external field, which couples to the order parameter. Indeed, in the linear 
response regime QCPs are typically characterized by singular susceptibilities at low frequencies [J, One can expect 
that the universality persists even beyond the linear response regime as long as the systems remains close to the 
criticality. Another possible situation where one can expect universal behavior is when the system is slowly driven 
through a QCP. In this case, since the dynamics is nearly adiabatic, one expects that low energy excitations will 
play the dominant role. Moreover, one can generally expect that non-adiabatic effects will be especially strong near 
singularities like QPTs, so that the dynamics will be dominated by the universal regime. This is indeed the case at 
least in sufficiently small dimensions [^, [T^ • The third situation where universality of dynamics can be expected 
is the response to fast small amplitude quenches (sudden changes in the tuning parameter) near the critical point. 
Analysis of sudden and nearly adiabatic dynamics near QCPs will be the central subject of this chapter. 
^ Typically analysis of slow, nearly adiabatic dynamics, is complicated by the fact that the usual perturbative 
approaches fail. It has been recently realized that adiabatic perturbation theory can become an efficient tool for 
analyzing the behavior of various thermodynamic quantities like the density of quasi-particles and the non-adiabatic 
energy (heat) generated during the process Although to leading order adiabatic perturbation theory often fails 

to give accurate quantitative analysis of various observables, it does give their correct scaling behavior with the rate of 
change of the external parameter. The advantage of this method is that it allows one to effectively reduce a dynamical 
problem to the static one and use the developed machinery for equilibrium quantum phase transitions. In this chapter 
we will demonstrate how this approach reproduces the correct scaling behavior of the response of the system both 
to sudden quenches of small amplitude and to slow nearly adiabatic sweeps through the critical point. We will also 
discuss the close connection between universal scaling of the density of quasi-particles and the scaling behavior of 
the fidelity susceptibility near the QCP, the quantity characterizing the overlap of the ground state wave functions 
corresponding to different coupling constants. In particular, we will show that the Kibble-Zurek scaling [3, [l^ can 
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be understood using this concept. We will illustrate some of the results using the transverse-field Ising model. In the 
end we will briefly discuss connections between adiabatic dynamics near critical points and dynamics associated with 
space-time singularities in the metrics, which naturally emerge in such areas as cosmology and string theory. 

I. BRIEF REVIEW OF THE SCALING THEORY FOR SECOND ORDER PHASE TRANSITIONS 

Continuous quantum phase transitions in many respects are similar to classical second order phase transitions. 
The main difference is that the quantum transition from one phase to another is driven by quantum rather than by 
thermal fluctuations, arising from the zero-point motion. So QPTs can happen at zero temperature. Examples of 
models where QPTs take place include quantum Ising and rotor models, the sine-Gordon model, various transitions 
from glassy to ordered phases in disordered systems and many others [l| . Many examples of QPTs will be discussed 
in consequent chapters in this book. Very often QPTs in d-dimensional systems can be mapped to thermal (classical) 
transitions in d + z dimensions, where z is the dynamical critical exponent One of the most important properties 
of both quantum and classical second order transitions is the universality of the low-energy long-distance properties of 
the system. This universality implies that the details of underlying microscopic models are not important near critical 
points. Instead the properties of the system can be well characterized by the parameter describing its proximity to 
the critical point (tuning parameter) and by universal critical exponents describing the singular behavior of various 
quantities with this parameter. We note that even though QPTs strictly speaking occur only at zero temperatures, 
the universality of the scaling governed by QCPs extends well into the finite temperature domain Recently there 
has been considerable interest in unconventional phase transitions, the description of which requires deviations from 
the standard framework. Some of these unconventional transitions will be discussed in other chapters of this book. 
Here we will not consider them since their dynamics is not yet well understood. 

A key quantity characterizing continuous phase transitions is the correlation length ^(A), which defines the length 
scale separating the qualitatively different behavior of e.g. spatial correlation functions of the fluctuations of the order 
parameter. This length scale diverges with the tuning parameter A as 



where v is a critical exponent. Divergence of the length scale is accompanied by divergence of associated time scales. 
For classical phase transitions the corresponding time scales are associated with relaxational dynamics. In quantum 
systems by the uncertainty principle a divergent time scale is also characterized by a vanishing energy scale A(A). 
The latter is usually associated with a crossover in the quasi-particle spectrum. E.g. it can represent a gap, or an 
energy scale where the dynamical exponent changes. Scaling of A defines another dynamical critical exponent z: 



There are many other critical exponents, however, these two will play the most important role in our discussion. To- 
gether with critical exponents, one can introduce scaling functions which describe long-distance low-energy properties 
of the system. For example, if we are talking about a QPT in a spin system then near the critical point we expect 
the equal-time correlation between the spin s at different positions to scale as 



where a is the scaling dimension of the spin s and F is some scaling function which approaches a constant when 
x/^ ^ 1. (Note that the divergence of the correlation functions at very small x will be cut off by non- universal 
short-distance physics). 

Other important quantities characterizing continuous phase transitions are susceptibilities, which describe the 
response of the system to external perturbations. For example, for a spin system the magnetic susceptibility describes 
the response of the magnetization to a small modulation of the magnetic field. From standard perturbation theory 
it is well known that susceptibilities are closely related to correlation functions Near QCPs static susceptibilities 
usually have singular non-analytic behavior characterized by their own critical exponents (see e.g. Ref. 0). Recently, 
it was realized that a very useful measure to analyze quantum phase transitions is fidelity susceptibility (FS) X/ or 
more generally a quantum geometric tensor (l6j . As we will see, this equilibrium concept and its generalizations will 
be very important for us later when we analyze dynamics near QCPs. Formally in a d-dimensional system FS is 
defined as 



C(A)^1/|A-A,r, 



(1) 



A(A)^1/^(A)^^|A-A,| 
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where |vE'„(A)) denote instantaneous eigenstates of the Hamiltonian H{X) and En{X) are the instantaneous energies. In 
this chapter we are only concerned with a non-degenerate ground state. FS appears in the leading order of expansion 
of the overlap of the ground state functions |(vl'o('^ + i5A)|\E'o(A)p in the powers of i5A: 

|(vI/o(A + a)|*o(A)P « 1 - SX^L^XfiX). (5) 

In the case when the system is translationally invariant and the operator d\H is local, i.e., d\H = J d'^xv{x), where 
X is the discrete or continuous coordinate, Eq. (j4]) implies that the scaling dimension of Xf is dim[xy] = 2Ay — 22 — d, 
where Ay is the scaling dimension of v{x) (i.e. at the critical point (v(x)v(x')) ~ _J_/|a; — x'P^"). This scaling dimension 
implies that the singular part of the FS near the critical point behaves as [la, [l3| 

X/(A)~|A-A,r(2A„-2.-rf), (6) 

In finite size systems the singularity of x / at the critical point is cutoff by the system size so that 

X/(A,)^L-2A"+2^+'^. (7) 

The crossover between the scalings ([6]) and (O occurs when the correlation length ^ becomes comparable to the 
system size L. For relevant or marginal perturbations the product Xv{x) should have the scaling dimension of the 
energy density: d + z. Because the scaling dimension of A is by definition 1/v we find that in this case A^, = 
d + z — l/f [3, Then Eqs. ^ and ([7]) for the singular part of Xf reduce to 

On top of the singular part there can be a non-singular part, which is generally non-sensitive to the proximity to the 
critical point. From Eq. (jSJ we see that the singular part of Xf gives is dominant near the QCP when diy < 2. The 
fidelity susceptibility can be expressed through the imaginary time connected correlation function p^] : 

X/(A) = ;^y^ dTT{{d,H{T)d,H{0))-{d,H(Q)r), (9) 

where d\H{T) = exp[_ffT]9A_ff exp[— _ffr]. The above relation immediately follows from the Lehmann's representation. 
From Eq. (jH) it is obvious that the singular scaling of Xf near the QCP is related to the power-law long-distance 
asymptotics of the correlation function. 

For the reasons which will become clear later is very convenient to generalize the fidelity susceptibility (|3]) to the 
generalized adiabatic susceptibility of order m [1^ : 

For m — 2 we recover the fidelity susceptibility: X2 = Xf- For m = 1 the corresponding susceptibility describes the 
second order correction to the ground state energy. The scaling dimension of Xm immediately follows from that of Xf- 
dim[xm] = dim[x/] — z{m— 2). The generalized adiabatic susceptibility can also be expressed through the connected 
correlation function: 

XrnW - drr"^-' {{d,H{r)d,H{Q)) - (9,i?(0))2) . (11) 

Example: transverse-field Ising model. Let us illustrate some of the generic properties mentioned above using the 
transverse-field Ising model, which is perhaps one of the simplest models showing quantum critical behavior. The 
Hamiltonian describing this model is 

Hi = -JY.r{9^f^(^M+i), (12) 

where trf and trf are the Pauli matrices, commuting on different sites, J is the energy scale for spin-spin coupling, and 
Jg is the energy scale for the transverse field. This model undergoes two QPTs at g = ±1. At small magnitude of 
the transverse field g the spins are predominantly magnetized along the z-axis, while at large \g\ they are magnetized 
along the direction of the magnetic field. In both cases there is a finite gap to the lowest energy excitations which 
vanishes at the critical point !1|]. This model can be solved exactly using the Jordan- Wigner transformation to map 
it to noninteracting fermions. 
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The Hamiltonian is diagonal in momentuni space and can be rewritten as 

Hi = j2^k{ihk-m, (13) 

k 

where ek — 2y/l + g'^ — 2g cos k and 7fc,7|. are the fermionic operators. The ground state of this Hamihonian, which 
is the vacuum of 7^, is written as a direct product: 

1*0) - (g) [cos((?fc/2) + isin(0fe/2)4cLj |0), (14) 

fc>0 

where |0) is the vacuum of original Jordan- Wigner fermions and ta.ndk = sin(fc)/(cos(fc) — g). Further details can be 
found in Sachdev's book [l|. 

It is convenient to introduce a tuning parameter A = g — 1 so that the transition occurs at A = Ac = 0. Near 
the critical point the expressions for 0k and dispersion simplify: tanOk sa —k/{k'^ + A) ss — fc/A and tk ~ 2-\/ A^ + k"^ . 
This scaling of energy immediately suggests that the critical exponents here are z — v — \. Indeed, according to 
our general discussion (see Eq. ^) the characteristic energy scale, which is obviously the gap in our case, scales as 
A = lAI^*^, which implies that zv = 1. Also, the spectrum clearly has a crossover from a constant to a linear function 
of momentum at k* ~ |A|, suggesting that there is a characteristic correlation length scales as ^ ~ 1/k* ~ 1/|A|; this 
implies that v ^ 1. Using the factorization property ([Ti)) it is easy to check that 

fc>0 ^ ^ fe>0 ^ 

In the thermodynamic limit for A 3> 1/i we can substitute the sum with the integral and find 

( \ \ f dk k \ / 1 

At the QCP and finite L Eq. (fT5| gives x/(0) ~ Both asymptotics are consistent with the expected scaling ([8|). 



II. SCALING ANALYSIS FOR DYNAMICS NEAR QUANTUM CRITICAL POINTS 

Above we gave a brief overview of the universal aspects of the equilibrium ground state properties of quantum 
critical systems. Our next goal is to see how similar universality emerges in out-of-equilibrium situations where the 
parameter A is tuned in time through a QCP. We expect that the dynamics will be universal if it is dominated by 
the low energy excitations generated near the critical point. In this chapter we will consider two specific situations, 
where such universality is found: (a) instantaneous quench of the small amplitude starting precisely at the critical 
point; and (b) slow quenches, where the tuning parameter changes gradually across the QCP. We will also allow for 
slow quenches to start (end) exactly at the QCP. These two situations correspond to two very generic scenarios. The 
first one is realized when the rate of change of the tuning parameter is fast compared to other relevant time scales. 
Though starting exactly at the critical point can require fine tuning, the scaling results will remain valid as long as the 
system is sufficiently close to the critical point both before and after the quench. The second situation is applicable 
to regimes in which the tuning parameter changes slowly in time. As we will show below these two scenarios can be 
actually described within the same unifying framework. In this section we will use qualitative arguments based on 
the scaling of the FS near the critical point. In the next section we will derive these scaling relations more accurately 
using adiabatic perturbation theory. 

We first start from instantaneous quenches. As we mentioned above, we assume the system is initially prepared in 
the ground state of some Hamiltonian Hq corresponding precisely to the QCP. We then suddenly apply perturbation 
XV, where A is a small parameter and V is some operator independent of A. We will assume that y is a relevant 
perturbation which drives the system away from the QCP. Let us discuss the expected scaling of various quantities 
using ordinary perturbation theory. According to general rules of quantum mechanics, the ground state wave function 
after the quench will be projected onto the basis of the new (quenched) Hamiltonian. Within ordinary perturbation 
theory the amplitude Sipn to occupy the excited state \n) is 



(17) 
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where the matrix elements and the energies are calculated for the unperturbed Hamiltonian. It is clear that quantities 
which commute with the new Hamiltonian like the excitation probability or the energy should scale as a square power 
of A. This follows from the fact that such quantities must be quadratic in Sipn- For example, the probability of 
exciting the system is 

P,, = ^ IS^^l^ = ^ = A^i^X/(0), (18) 

where the prime over the sum implies that the ground state is excluded from the summation. The excitation probability 
is meaningful only for quenches with the amplitude vanishing in the thermodynamic limit such that ^(A) > L). For 
larger amplitude quenches Pex ~ 1. In these situations one needs other quantities to characterize the system's response 
to the quench. One of such quantities can be diagonal entropy [13]: Sd — — X] IV'nP In iV'nPj where IV'nP are the 
probabilities to occupy eigenstates of the quenched Hamiltonian. Clearly Sd is a possible measure of the non-adiabatic 
response of the system, which is well defined for quenches of all amplitudes. Another possible measure characterizing 
the non-adiabatic response of the system is the energy, or more precisely the non-adiabatic part of the energy change 



due to transitions. This excess energy, or simply heat 21|, is the difference between the energy after the quench and 



the instantaneous ground state energy. Within the perturbation theory we find that the heat density Q is 

Q(^) = ^d E(^" - ^o)|5^.P = T^A^ E ^ A\MO), (19) 



where we introduced the energy susceptibility X-e(A) = Xi(A)- If we additionally assume that excitations in the 
system are carried by well-defined quasi-particles, i.e., that the energy eigenstates |n) are also the eigenstates of the 
quasi-particle number operator (which is usually the case only in integrable models), then we can also compute the 
density of quasi-particles, 

ncx(A) = ^^E^"!*^^"!' = (20) 

n 

where A^,^ is the total occupation of each eigenstate n 

Xn-^E^"(^— (21) 

In most systems the number of quasi-particles is actually not conserved due to various collision processes. However, 
the quasi-particle number can be still a very useful quantity if we are interested in times shorter than the relaxation 
time, i.e., times in which collisions do not play any significant role. We see that in all these situations we expect 
^ox(A), (9(A), and ncx(A) to scale quadratically in A to lowest order. The scaling for the d-entropy can ac quir e 
additional logarithmic corrections. We will not discuss it here for brevity since it is very similar to that of tIcx il9|. 
We note that this quadratic scaling should be contrasted with the usual linear response relations where various 
quantities are proportional to the first power of A, e.g., 

(V'(A)|T/|V'(A)) w const + 2Ay^^^^ffl^ ^ const + 2i'^Axis. (22) 

To obtain linear dependence it is obviously crucial that the observable of interest does not commute with the Hamil- 
tonian in the excited state. Such situations can also be analyzed. However, we will not consider them. 

It is very important to remark that we are dealing with extended systems. As we saw in the simple example of 
the transverse-field Ising model the FS is extensive in the system size ()16p . This situation is generic if we consider 
global (spatially-uniform) perturbations. Strictly speaking this implies that the validity of the perturbation theory is 
restricted to very small perturbations, where Pcx(A) <C 1. However, we know very well that perturbative approaches 
typically have much larger domains of applicability. It is usually important that only changes of intensive quantities 
like the energy per unit volume or the density of quasi-particles remain small. Indeed, the probability to produce the 
next quasi-particle excitation is not affected much by the presence of other quasi-particles if those are very dilute. Thus 
we expect that Eqs. p9|) and ((20|) have a much broader domain of applicability than Eq. ([T8|) (see also Refs. [H, [l^ 
for additional discussion). 

As we discussed earlier, near QCPs susceptibilities may become divergent (see e.g. Eq. (IT51) for the Ising model) or 
acquire some other type of singularity. Such divergences can invalidate the quadratic perturbative scaling of various 
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observables and these situations will be the next point of our discussion. Let us look into e.g. Eq. (|18p for the 
probability of exciting the system. Near the critical point X/('^) may diverge as A — >■ with the divergence being 
cutoff by the system size (see Eq.. [8]). This results in the superextensive scaling of the transition probability: 

Pex(A) ^ X^L^'" (23) 

for 2/v > d. It is clear that Pox should saturate at unity when A becomes bigger than A* ^ l/L^/"^, which is 
accidentally equivalent to the condition ^(A) < L where the fidelity susceptibility changes its form (see Eq. [S]). 
Likewise for 2/ v > d + z for the heat density we find 

Q(A) - x^L^/'^-d-- (24) 

for A < A*. However, unlike the probability of excitations, we can extend the perturbative expression for the heat (jl9p 
to the regime A > A* by simply evaluating the susceptibility xe a-t the final value of the quench parameter. More 
accurate derivation based on the adiabatic perturbation theory, which we discuss in the next section, shows that one 
needs to change xb(0) ~^ Jo XE{X')dX' in Eq. (fT9|) . However, the scaling of this integral is the same as the scaling of 
Xe{,X). Then we find that for A > A*: 

Q(A) - \X\^'^+'>. (25) 

Finally if we additionally assume that dominant excitations are coming from isolated quasi-particles (which is often 
the case) we find that for A > A* 

nex - |A|'^^ (26) 

while for A < A* the scaling for riox is similar to that of Pcx given by Eq. 

Of course, we expect these scalings to remain valid as long as the corresponding exponents do not exceed two. 
Otherwise, the low-energy singularities associated with the critical point become sub-leading (corresponding suscep- 
tibilities do not diverge) and the perturbative quadratic scaling is restored (though the singularities still appear in 
higher order derivatives with respect to A). 

There is a very simple argument, which reproduces the scalings (1251) and (|26p . A quench of amplitude A gives us 
natural length and energy scales ^ ^ 1/1 A]"^ and A ~ jAj^". We thus might expect that the quasi-particle excitations 
with energies larger than A, or equivalently with momenta larger than 1/^, will not be much affected by the quench 
and can be treated perturbatively, yielding the quadratic scaling. At the same time, for states with energies less 
than A the quench will be effectively very strong, so that they will be excited with the probability of the order of 
unity. Thus the density of quasi-particles will scale as l/C' ~ X'^'^ and we reproduce the scaling ((26l) . li dv <2 then 
this contribution to riox will be dominant over the perturbative contribution coming from high energies (because it is 
proportional to smaller power of A). Likewise, we can obtain the scaling of heat (|25p by noting that the low energy 
excitations carry energy of the order of A. Thus their contribution is Q ~ noxA ~ jAj'-''"'"^''''. When the exponent 
(d -I- z)v becomes more than two this low energy contribution becomes sub-leading and the perturbative quadratic 
scaling is restored to leading order in A (as we mentioned singular non-analytic terms can still persist to higher order 
in A). There is a related and very intuitive way of deriving the scaling (j26l) adopting the arguments of Kibble and 
Zurek ^^'^ situation of sudden quenches. Namely, we can interpret the scale ^ as a typical distance between 

generated quasi-particles. Then the quasi-particle density is ~ \X\'^^ . 

It is interesting that Eq. (|23p can also be understood purely by symmetry arguments. Let us assume that the 
phase corresponding to finite A is characterized by some broken symmetry. The symmetry becomes well defined when 
^(A) ~ 1/|A|'' becomes less or comparable to the system size: ^(A) < L. This defines the minimal quench amplitude 
A* ^ l/L^/" at which the symmetry is formed. Since the critical point does not correspond to any broken symmetry 
we anticipate that the overlap between ground state wave-functions \'4>{0)) and |'0(A)) vanishes for A > A*. This 
implies that 

Xfm\X*f > 1 (27) 

or equivalently x/(0) > L^/'^"'*, which is actually the correct scaling. In the case, dv > 2, Eq. ([77|) is satisfied even 
if xf(A) approaches a constant independent of the system size at the critical point. Note that the argument, in 
principle, allows X/(0) to vanish at the QCP for dv > 2 in the thermodynamic limit. However, as we discussed above, 
generally one can anticipate that x/(0) is nonzero due to high-energy non- universal contributions not sensitive to the 
presence of the critical point. 

Let us now consider a somewhat different setup, where we still start at the critical point but instead of suddenly 
quenching the parameter A we assume that it gradually increases over time as a power law A(t) = vt^ /r\Q(t), where 
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the factor of l/r! is added for convenience. Here Q{t) is the step function and i; is a real constant, which controls 
the proximity to the adiabatic limit. For linear quenches r = 1 the parameter v plays the velocity of the quench, for 
quadratic quenches r = 2 it is the acceleration and so on. It is interesting that in the limit r — )■ this parameter v 
becomes the quench amplitude. For simplicity we also assume that for r > the final value of A is sufficiently far 
from the critical point. As in the case of sudden quenches, it is instructive first to perform perturbative in v analysis 
of the new state. Using adiabatic perturbation theory Oll^], which we will briefly discuss in the next section, one 
can show that to leading order in v the probability to occupy the state \n) in the instantaneous (co-moving basis) is 

I'^^"! ^ (i?„-£;o)2-+2' 

where as before all matrix elements and energies are evaluated at the critical point. We note that the same expression 
applies to the opposite situation, where one starts far from the critical point and changes the coupling linearly in time 
until the critical point is reached. This expression can be also generalized to situations when both initial and final 
couplings are finite [Tsl] . Then the total probability to excite the system reads 

^-(-) - E Jt^l^^^ - v^L'^X2rM0), (29) 

where X2r+2(A) is the generalized adiabatic susceptibility introduced in the previous section ()10|) . One can similarly 
introduce susceptibilities describing the heat density Q and the quasiparticle density n^x for slow quenches. For 
example, it is easy to see that if the quench ends in the critical point then 

Q(w)«w'x2,.+i(0). (30) 

If, on the other hand, the quench ends far from the critical point then the energy of the system and the number of 
quasi-particles are evaluated in the states adiabatically connected to the states |n) at the critical points, thus the 
expressions for Q and n^x become in general nonuniversal. The situation simplifies in integrable models, where e.g. 
the quasi-particle number is conserved during the adiabatic evolution. The same is true in nonintegrable models if 
the excitations are topological and have very slow relaxation. The ambiguities of choosing the basis also disappear if 
we consider a cyclic process where coupling first linearly increases in time, then saturates, and then decreases back 
towards its original value, e.g. X{t) oc t^{tf — tYQ{t{tf — t)). In this case, the corresponding susceptibilities essentially 
will have an extra factor of two. Also, there are no complications with choosing the basis if we consider a reverse 
process where we stop at the critical point. 

Now we are in the position of more closely analyzing the scaling of Pox and other quantities. Like in the case of 
instantaneous quenches X2r+2(-^) for r > can diverge at the critical point because of low energy contributions to 
Eq. (1^^ . Since as we discussed earlier dimx2r+2 = d—2/i'—2zr we see that this is indeed the case for d < Ijv + lzr. 
As in the case of sudden quenches this divergence leads to super-extensive scaling for the excitation probability: 

P,x(v) ^ |t;pi2/-+2-. (31) 

This probability becomes close to one at the crossover rate v* ~ l/L^/''+^''. For faster quenches, v > v* ^ the excitation 
probability is no longer informative and we need to consider other quantities. If we assume that the quasi-particles are 
well defined then in the case \ v\ > \v*\ we can extend Eq. ([^ as we did for sudden quenches, namely evaluating X2r+2 
at some characteristic coupling slightly away from the critical point. Noting also that = dim[A] = dim[t;] — zr we 
see that the scaling dimension of v is dimfu] = zr + l/v. This means that the characteristic value of A, which should 
enter the susceptibility is related to v via zvr). Therefore we recover that for \v\ > \v*\ 

nox~|w|^, (32) 

which for linear quenches r = 1 is indeed the correct scaling first suggested in Refs. Q and later generalized to 
nonlinear quenches in Refs. [l^ll^]. If we assume that in the final state the spectrum is gapless and characterized by 
the exponent z (e.g. if the final state of the evolution corresponds to the critical point), we find the the heat is also 
universal. In particular for \v\ > \v*\: 

Q^\v\^. (33) 

Equation (|33p is also the correct scaling first suggested in Ref. for r = 1. As in the case of sudden quenches 
these scaling results are expected to be valid only if the corresponding exponents in the powers of v are smaller than 
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two since otherwise the corresponding susceptibihties X2r+2(A) or X2r+i(A) are not divergent genericaUy having only 
cusp singularities at the QCP. In this case the corresponding non-analytic asymptotics of Pcxi and tt-cx become 
subleading. 

There is also a very intuitive explanation of the scaling p2l) . As we mentioned in the previous section, critical 
points are characterized by the quasi-particle energy scale A ^ jAj^"^ — \vtY^ (we note that the scale A, relevant 
to our discussion, is always associated with quasi-particle excitations, even if quasi-particles are ill defined; many- 
body energy levels are generally exponentially close to each other). If this energy scale changes sufficiently slowly 
in time then the energy levels have time to adjust to this change and adiabatically evolve. However, if A changes 
sufficiently fast then the adiabaticity breaks down and the states are excited. To find the crossover energy scale 
separating adiabatic and diabatic states we can use the simple Landau-Zener-Majorana-Stiickelberg (LZMS) criterion 
dK/dt ^ A^. Using two facts that A oc \X\^'' and A ~ vV we find that A ~ \v\'''' / i^^^+i) _ r^j^^^ ^^^^j^ corresponds 
to the characteristic momentum k ~ and the characteristic length scale ^ ~ 1/fc. The number of excited 

quasi-particle states in the spatially uniform system is then riox ^ fc"^ ~ \^y]d.u / {zvr+i) ^ which is exactly as in Eq. (|32p . 
We note that in this form the argument does not require that the initial or final point of the evolution coincides 
with the critical point. It is only important that the QCP is crossed during the time evolution. In the next section, 
in which we discuss adiabatic perturbation theory, we will show that this is indeed the case. The crossover to the 
quadratic scaling is also expected to be generic. In general in the expressions and (15(11) one has to evaluate the 
corresponding susceptibilities in the initial and the final points of the evolution. Away from the phase transitions 
we expect no singularities in the corresponding susceptibilities and thus the quadratic scaling will hold. This simple 
scaling argument can be also reformulated in the spirit of Kibble- Zurek mechanism [l3.[l5j. Namely, one can interpret 
the length scale ^ as a characteristic distance between generated quasi-particles and the inverse energy scale A~^ as 
a time scale where the process is not adiabatic and the defects are generated. 

Example: transverse-field Ising model and multi- dimensional extensions. Let us now illustrate how the scalings we 
derived apply to the specific example we introduced in Sec. H] We start from sudden quenches. We also focus on 
the thermodynamic limit i — > cx) and A is independent of L. Using the explicit structure of the ground state wave 
function (IT4)) it is easy to check that the overlap of two different ground states corresponding to a particular pair of 
fermions with momenta fc, —k is 



(vI/g(A)|*?;(A'))=cosp,-0l.)/2], (34) 

where 9k and 0'^, correspond to the couplings A and A' respectively. This implies that the probability of exciting a pair 
of quasi-particles with momenta k and —k by quenching the parameter A' to A is Pox(fc) = sin^ {[9k — 0^]/2). Noting 
that we are interested in the limit A' = and using low energy expressions for Ok and O'k we find 



1 - 



Vfc2 + A2 



(35) 



The density of quasi-particles excited in the quench is then obviously obtained by integrating Pcx{k) over different 
momenta 



(36) 



This scaling indeed agrees with our general expectation Ucx ^ |A|'''', noting that d = ly = 1. Similarly we can find 
heat, noting that each quasi-particle with the momentum k carries energy « 

2VA2 + fc2. Then 



r dk r- 



1,9, 27r 



2^AVx(fc)«^A2lnj^, (37) 



which also agrees with Q ^ | Aj^''"'"^^'^. However, because in this case the power [d-\- z)v is exactly equal to two, we find 
an additional logarithmic dependence on both A and the cutoff tt. The logarithmic dependence is natural at the point 
where we expect crossover from the exponent two at [d -f z)v > 2 to the exponent less than two in the opposite case. 
Even though the transverse-field Ising model is defined only in ID, one can extend it to higher dimensions by formally 
considering the free-fermion Hamiltonian (|13|) in higher dimensional lattices. E.g. in 2D a very similar Hamiltonian 
represents the fermionic sector of the Kitaev model [1^ or describes electrons and holes in graphene [l^. Then it is 
easy to check that rtox has quadratic scaling with A in 2D and above with extra logarithmic corrections in 2D where 
dv = 2. The heat has quadratic scaling above ID (see Ref. [l^ for additional details). 

The transverse-field Ising model can be also solved for slow quenches. For simplicity we will analyze only the linear 
dependence A = vt. Note that the wave function factorizes into a direct product of states corresponding to different 
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momenta with either zero fermions or two fermions in each state. Thus the dynamical problem factorizes into a direct 
sum of LZMS problems. If the magnetic field linearly crosses the QCP then the transition probability is given by [23] 

PoAk) ~ exp(-27rfcVw) (38) 

The density of quasi-particles generated in such process is obtained by integrating Pcx(^) over different momentum 
states yielding in the slow limit riox « ^/v / {2t: \/2) . This is indeed the expected scaling n^^ ~ |^;|<i'^/(2'^+i) for a 
particular set of exponents d = v = z = I. The problem can be also solved for the linear quench starting at the 
critical point giving identical scaling with a slightly different prefactor [l^. To see the crossover to the quadratic 
scaling one needs to consider extension of the Hamiltonian (|13p to higher dimensions. In this case for each momentum 
state we are effectively dealing with half LZMS problems, i.e. LZMS problems where the initial coupling corresponds 
to the minimum gap. Asymptotically, in the slow limit f <C e^, this probability is found from Eq. (|28p : 

Pex(fc) « (1024)-it;Vfc* (39) 

It is clear that above four dimensions the density of quasi-particlcs ^ J d'^kp^x will be quadratic, dominated 
by excitations to high energies of the order of the cutoff. Below four dimensions the integral converges at large k 
so transition probabilities are dominated by small momenta k ^ ^ and the scaling ricx ^ — \v\^^'^ is 

restored. 



III. ADIABATIC PERTURBATION THEORY 
A. Sketch of the Derivation 

In this section we will present a more accurate derivation of the scaling of various quantities near the critical 
point. In particular, our aim is to derive the expression for the transition probability in the leading order in v. 
A very convenient framework to analyze these scaling laws is given by adiabatic perturbation theory. As we will 
see this approach will allow us not only to reproduce Eq. but also extend this result to the regime where the 

corresponding susceptibility X2r+2(A) diverges at the QCP and justify Eqs. (|32p and (|33p . In Ref. [s*] this theory 
was originally applied to derive the scaling of the density of quasi-particles p2p . In Ref. [llj it was shown that 
this approach correctly predicts the crossover between analytic and non-analytic regimes of scaling of heat (excess 
energy) with the rate v for linear quenches. In Refs. [l^ [l9| this theory was extended to sudden quenches near 
critical points and also reproduced correct scaling results. While adiabatic perturbation theory is not quantitatively 
accurate in non-analytic regimes (where the response is not quadratic) i-e., it does not correctly reproduce 

the prefactor, it predicts correct scaling laws in many different situations. The only known exceptions are in the 
so-called non-adiabatic regime, where the system size or other macroscopic length scale enters the scaling of heat or 
quasi-particle density Such regimes can appear e.g. if we are dealing with low-dimensional systems, which 

have low energy bosonic excitations, especially at finite temperatures, where the violation of scaling p2p comes from 
the overpopulation of low-energy modes. This regime is quite special and we will not consider it here (we refer the 
reader interested in more details to Refs. [l9|). 

In the beginning of the section we will closely follow the discussion of Refs. 0, [2^. We consider a very general 
setup where the system is described by a Hamiltonian Hit) = i/o + '^(i)^j where A(t) monotonically changes in time 
between initial and final values Ai and A/. For simplicity in this chapter we will focus only in linear dependence 
A(t) = ut, though this assumption is not important. The limit u — )■ cx) corresponds to a sudden quench and f — > 
to a slow quench. We will always assume that for sudden quenches A/ is close to A^, while for slow quenches this 
condition is not necessary. In both cases adiabatic perturbation theory will be justified by the proximity of the system 
to the ground state after the quench. 

Our goal is to approximately solve the Schrodinger equation 

%dt\i=)=H{t)\^), (40) 

where is the wave function. It is convenient to rewrite Eq. (j40p in the adiabatic (instantaneous) basis: 

\m)^Y.n^n{t)\n{t)), H{t)\n{t))=Er.{t)\n{t)), (41) 

where £'„(t) are the instantaneous eigenvalues. The eigenstates \n{t)) implicitly depend on time through the coupling 
X{t). Substituting this expansion into the Schrodinger equation and multiplying it by (m| (to simplify our notations 
we drop the time label t in \n(t))) we find: 



(42) 
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Next we will perform a unitary transformation: 

an{t) = an{t) exp [-iQn{t)] , Quit) - Jl En{T)dT. (43) 

The lower limit of integration in the expression for Qnit) is arbitrary. We chose it to be equal to ti for convenience. 
Then the Schrodinger equation becomes 

an{t) = - E™ a^{t){n\dt\m) exp [i(e„(i) - Q,n{t))] ■ (44) 
In turn this equation can be rewritten as an integral equation 

an{t) = dt'E^a™(i')(n|at'|m)e^(«^"(*')-«"(*')). (45) 
If the energy levels En{t) and E„i{t) arc not degenerate, the matrix element {n\dt\m) can be written as 

{ \d \ ) {nl^^tHlm) ^^^^ {n\V\m) ^^^^ 

En (t) — Em (t) En (t) — Em (t) 

If X{t) is a monotonic function of time, like in our case, then in Eq. (I45p one can change variables from t to X{t) and 
derive 

«n(A) =-lt rfA'E„a™(A')(n|5v|™)e^(«"(^')-®-(^'», (47) 

where 

QniX)^ J^^dX'EniX'KX')-\ (48) 

Equations (P5|) and (|T7)) suggest a systematic expansion in the transition amplitudes to the excited states. We are 
interested in either the limit A = u — ?► 0, which suppresses transitions because of the highly oscillating phase factor, 
or the limit of small — Ai|, where the transitions are suppressed by the smallness of the integration domain. To 
leading order in adiabatic perturbation theory only the diagonal terms with m = n should be retained in the sums in 
Eqs. (pSj) and (|T7)) . These terms result in the emergence of a Berry phase (29j : 

^n{t) - -tJl dt'{n\dt,\n) = -z/^';^*^ dX'{n\dx,\n) , (49) 

so that a„(i) ~ a„(0) exp[— i<I>„(i)]. In general, the Berry phase can be incorporated into our formalism by performing 
a unitary transformation a„(t) a„(t) exp[— i$„(t)] and changing 0„ — )■ 9„ + $„ in Eqs. (1451) and (|47|) . 

In many situations, when we deal with real Hamiltonians, the Berry phase is identically equal to zero. However, in 
some cases when more than one coupling constant change in time, the contribution of the geometric phases can be 
important, so that it can change the results for scaling of the physical quantities near the phase transition (s^l- Here 
we only note that since the geometric phase is related to the topology of the phase space, the evolution of physical 
quantities can depend on the path in the parameter space. The effects of the geometric phase are enhanced near 
the diabolic points corresponding to the level crossings. Geometric phase effects can be important for open systems, 
which effectively can be modeled by non-Hermitean (complex) Hamiltonians. 

Assuming that the geometric phase is not important, let us compute the first order correction to the wave function 
assuming that initially the system is in the pure state n = 0, so that ao(0) = 1 and a„(0) = for n ^ 0. To leading 
order in A we can keep only one term with to = in the sums in Eqs. (1451) and (I47p and derive 



anit) di'(n|9t,|0)e'[«"(*')-«5o(t')]. (50) 

or alternatively 

Q„(A) « -/;^_ dA'(n|9A'|0)e*[®"(^')-®°(^')l. (51) 

The transition probability from the level |0) to the level \n) as a result of the process is determined by |ara(A/)p. Let 
us note that in the limit of A = u — >■ one can expect two types of contributions to Eq. ([?T|) . (i) A non-analytic 
contribution comes from the saddle points of the phase difference Qn{X) — Qm{X), which in turn correspond to the 
complex roots of i?„(A) = Eq{X). These terms result in exponential dependence of the transition probability on 
rate jap ~ exp[— y4/|A|] like in the usual LZMS problem (see Ref. [l^ for the additional discussion), (ii) An analytic 
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contribution comes from the moments where we turn on and turn off the process. This second contribution, to leading 
order in A, can be obtained by integrating Eq. ([5T|) by parts [l3l[2^: 



a„(A/) 



iX 



{n\dx\0) 



E„i\)-EoiX) 



,j(e„(A)-eQ(A)) 



(52) 



Ai 



From this we find the analytic part of the transition probability 



\{n\dyM' 



\{n\d>.M' 



- 2v 



(£;„(A,)-So(A.))2 
2 {n\dxM 



{En{Xf) 

{n\dx,\Q) 



En{X^) - EoiX,) En{Xf) - Eo{Xf) 



cos [Ae„o 



(53) 



where A9„o = ©n 
will average out to 
off the coupling A. 
we recover Eq. 
evolution then we 
of X{t). It is strai 
expression can be 
perturbation theor 



(Xf) — 8o(A/) — 8„(Ai) + 8o(Ai). Usually if we deal with many levels the last fast-oscillating term 
zero. The remaining first two terms are the transition probabilities associated with turning on and 
If only the first term dominates the scaling (e.g. because the initial state corresponds to a QCP) 
I for r = 1 . If we deal with nonlinear quenches such that A = at the initial and final points of the 
need to continue integrating Eq. (1511) by parts until we hit the first non- vanishing time-derivative 
ghtforward to see that in this way we reproduce Eq. (j28l) for any integer power r. The resulting 
analytically continued to all positive r. Interestingly it correctly reproduces the result of ordinary 
y for sudden quenches r = where the parameter v plays the role of the quench amplitude. 



B. Applications to Dynamics near Critical Points 

Let us now perform the scaling analysis of Eq. ([51]) . First we consider the situation of sudden quenches with the 
initial coupling corresponding to the critical point A^ = (or alternatively A/ =0). We will assume that excitations 
are dominated by quasi-particles created in pairs with opposite momenta k,—k. As we saw in the previous section this 
assumption is well justified for non-interacting models like the transverse-field Ising model. In general, for interacting 
non-integrable models eigenstates consist of quasi-particles dressed by interactions and the quasi-particles acquire 
finite life time. If the quasi-particle nature does not qualitatively change because of the interactions we expect that 
the scaling of the quasi-particle density will remain the same as the scaling of the excitation probability.In fact, the 
whole scaling analysis can be performed in the spirit of Sec. [ill which does not rely on the assumption of having 
well-defined quasi-particles. However, such an assumption makes all derivations more transparent. 

Assuming that |n) corresponds to the quasi-particle pair with momenta k, —k (to simplify our notation we denote 
such pair states by |fc) and the corresponding energy of the pair by £k) we find 

Since the scaling dimension of the operator AV^ is the same as the scaling dimension of energy we expect the following 
scaling of the matrix element: 

<'l^'l»> -kwiAn, ,55, 



efc(A)-eo(A) A 



where G{x) is some scaling function. We anticipate that in the limit a; ^ 1 we have G{x) ~ x~^^'^ so that the matrix 
element becomes independent of A at fc 3> A"^. In the opposite limit, x 1, the scaling function can either saturate 
if there is a gap or vanish as some power of x if there is no gap. It is easy to check that for the transverse-field Ising 
model the scaling assumption is indeed satisfied with G{x) oc x/(x^ + 1). This scaling ansatz immediately allows us 
to analyze the behavior of the quasi-particle density: 

"ox~ -^EfclttfeP- (56) 

Then taking the thermodynamic limit where the sum over k becomes an integral and changing variables A — )■ Xfi] 
and fc — >■ |A/|''^ we immediately find 



^ex(A/)«|A/|'^'' 



V 



(57) 
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This expression gives the right scaling ([26)1 provided that the integral over ^ converges at large ^. The convergence 
clearly depends on the large f asymptotics of the function G'(^/|?7|'^) ~ l/j^l^/"^. The integral over ^ clearly converges 
provided that d <2/v or dv <2. In the opposite case the integral over momenta is dominated by large k » lA/l"^. 
In this case Eq. ([55]) reduces to the perturbative quadratic result, e.g. nox(A/) ~ 2AjX/(0), where the factor of two 
reflects that each excited state contributes to two quasi-particles. Note that the assumption of scale independence of 
the transition matrix element is equivalent to the assumption that Xf does not vanish at the critical point for dv > 2. 
As we pointed out earlier, this is generally expected because transitions to the high energy states are insensitive to the 
proximity to the critical point dominate excitations in the system (see detailed discussion in Ref. [l9| for a particular 
sine-Gordon model). Likewise, one can derive the correct scaling for the heat density (l25t provided that {d + z)!^ < 2; 
one reproduces the quadratic perturbative result in the opposite case. 

In a similar fashion we can derive the results for the adiabatic case, when the coupling constant changes slowly in 
time. For simplicity we focus again on linear quenches X{t) — vt. For now it will not be important whether we start 
(finish) at the critical point or the initial and final couplings are on different sides of the QCP. The expression for the 
transition amplitude then becomes 



akiXf) 



dX 



{k\V\0) 



efe(A)-eo(A) 



exp 



rfA'(efe(A')-eo(A')) 



(58) 



Let us now introduce another scaling function F{x) according to efe(A) — eo(A) — |A|^''F(fc/|A|'') with the asymptotic 
power law F(x) ~ \x\^ for \x\ 3> 1 (the small x limit depends on whether the system is gapless or not). A convenient 
change of variables is fc = \v\'^/^^''~^^^r] and A = jup/^^'^+^^r;. Then the expression for the density of quasi-particles 
becomes 



d'^k 
(2^ 



\v\ --+1 



d'^rj 



where 



a{ri) 



if 



ddc ( f 



exp 



(59) 



(60) 



An analogous expression is found for heat. If v is small then the limits of integration over f can be extended to 
{—CO, oo) if Aj < and A/ > or to (0, oo) if = and A/ is finite. Indeed, the integral over ^ is always convergent 
because of the highly oscillating exponent. If additionally the integral over the rescaled momentum rj can be extended 
to oo we get the desired scaling for the quasi-particle density (|32]) . 

To analyze the convergence (assuming for simplicity that |A/| 3> |Ai| and Xi is close to the QCP) we note that at 
large k the asymptotic expression for the transition amplitude a{k) according to Eq. ([5^ is 



a{k) 



IV 



{k\V\Q) V 



(61) 



where we kept only the term corresponding to the initial coupling A^ = because it describes the singularity. We 
thus immediately infer that the integral over the momentum k (and thus over the rescaled momentum 77) can be 
extended to 00 (so that the scaling ([55)) is valid) as long as d < 2zv + 2/v or equivalently dv/{zv -I- 1) < 2. In the 
opposite case, the usual adiabatic perturbation theory is restored and we recover the quadratic perturbative result 
^^ex ~ u^(x/(Ai) + X/(A/). We expect this scaling to be valid even if both Xi and A/ are away from the quantum 
critical region. Indeed, in this case one does not expect any divergences in susceptibilities Xm(A) at both A^ and 
A/ so the quadratic scaling should dominate the dynamics over the higher critical power emerging from 

the scaling argument. A somewhat special situation emerges if the relevant susceptibilities Xm(A) are very small at 
Xi and Xf. Then one might expect that the scaling |-ti|'*'^/(^'^+i) will be applicable even when the exponent exceeds 
two. A similar situation occurs if the initial and final couplings are finite but the time dependence of A(i) is smooth 
near ti and tf and linear only near the critical point. These expectations are well justified only for non-interacting 
systems like the transverse-field Ising model, where dynamics can be mapped to independent LZMS transitions. In 
general, interactions lead to dephasing of quasi-particles, which is equivalent to resetting the dynamical process each 
time the quasi-particle phase is lost. Then one can expect the quadratic scaling with v will be again restored for 
dv/{zv + 1) > 2. At the moment this issue remains an open problem. We note that in the case dv/{zv + 1) one can 
expect that the non-analytic term coming from low energies will still survive; it will just become sub-leading in the 
limit of small v. 

Another important point is that the susceptibilities X2r+2{X) for r > can diverge even away from QCP as long as 
the system is gapless. Indeed e.g. for bosonic theories gapless phases correspond to v ^ 00. Then the susceptibility 
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X2r+2 diverges for d < 2zr. This means that even in gapless noncritical phases in low dimensions the quadratic scahng 
can be violated and replaced by Hex ~ \v\'^/^^ and Q ^ |^;|(<i+^)/^''. This is consistent with general prediction of less 
adiabatic response in low dimensions [l^. It is also interesting to note that in the noncritical case for d/zr < 2 the 
density of quasiparticles (or for (d + z)j zr < 2 the heat density) depend only on the total time of the quench, but 
not its shape. Indeed since the total quench time r ~ we see that ricx ^ \v\'^^^^ ~ (1/t)'^/^, i.e. independent 

of r. On the other hand for d/zr > 2 and in gapped phases in arbitrary dimensions ticx ^ (1/''')^'^ is very sensitive 
to the exponent r and thus to the shape of the quench. The existense of these two regimes is true in other situations 
as well. Thus in Ref. [3l| the two regimes of sensitivity of the energy (heating) to the shape pulse were discovered 
for interacting gapless fermionic systems. In general relevance of the ramping protocol to the non-adiabatic dynamics 
in interacting many-body systems remains an open problem, which is likely very important for understanding the 
problem of dissipation. 



C. Quenches at Finite Temperatures, and the Role of Quasi-particle Statistics 

So far, we have exclusively focused on the situation where the system is prepared in the ground state. A natural 
question arises of how one can extend these results to the finite-temperature domain. It is well established that 
in equilibrium near a QPT influence of the QCP extends into finite temperatures in the so-called quantum critical 
region [I] which can be quite extensive. We can anticipate that a similar situation occurs for non-equilibrium systems. 
Let us point out that the definition of temperature can be somewhat ambiguous away from equilibrium. One can 
consider two natural setups where the system is coupled to some external reservoir even during a dynamical process 
(see e.g. Ref. 32]). In this situation results are sensitive to the coupling strength to the bath and in this sense are not 
universal. However, one can imagine a quite different setting, where the system is initially prepared in some thermal 
state and then the consequent dynamics is purely Hamiltonian, i.e., the effects of the bath are negligible during the 
time evolution. Such a setup was considered in Refs. [l^ [l9l . Issj . In this setup, no coupling to the environment is 
assumed during the evolution and one can expect universality of the results. Such setups routinely appear now in 
the cold-atom context in which systems, after initial preparation, are essentially isolated from the environment [s^l- 
They also appear generally in thermodynamics when one considers adiabatic or nearly-adiabatic processes (in which 
changes happen on time scales faster than the time of equilibration with the environment) (35j . Universality can only 
be expected in this setup if the system is initially prepared close to a QCP. So in this section we restrict our analysis 
only to such situations. 

On general grounds one can expect that at finite temperatures the statistical nature of low energy excitations 
should play a significant role. Indeed, if we are dealing with bosonic low energy quasi-particles we expect that 
the initial thermal population enhances the transition probability to the corresponding quasi-particle modes. The 
reason is that bosonic excitations have a bunching tendency, i.e., a tendency to stimulate transitions to a particular 
mode if it is already occupied. On the other hand, if quasi-particles are fermionic then Pauli blocking occurs due to 
already-occupied modes and fewer quasi-particles can be additionally excited. In principle, one can expect situations 
where low-energy excitations are either described by quasi-particles with fractional statistics like in Kitaev j25l] and 
sine-Gordon models [l^, or by many-particle excitations without well defined statistics. There the effect of finite 
temperature on dynamics is unknown at present and thus will not be discussed here. 

If quasi-particles in the quantum critical region are bosonic then it is straightforward to show that the number 
of excited quasi-particle pairs in the mode with momentum q at initial temperature T is related to the number of 
quasi-particle pairs created at zero temperature via a very simple expression valid for any dependence of X{t) fl2j: 

nex(<Z, T) = neq(g, T) + <(</) coth(eg/2r). (62) 

In Eq. (j62p . nc^{q^T) is the initial equilibrium population of the bosonic mode described by the Bose-Einstein distri- 
bution function, is the (initial) energy of the quasi-particle, and n^^(ci) is the number of quasi-particles created in 
the same dynamical process if the initial temperature is zero. At small temperatures 3> T this expression clearly 
reduces to the zero temperature limit. In the opposite high temperature limit the transition probability is enhanced 
by the factor 2T/e^ ^ 1. Similarly, in the fermionic case we have [l^ 

nox(g, T) = Ucqiq, T) + n"^{q) tanh(eg/2T). (63) 

This expression also reproduces the zero-temperature result for eq ^ T while yielding suppression of the transition 
probability by a factor Cq/ (2T) in the opposite limit. Since the total number of dynamically excited quasi-particles is 
found by summing Ucxiq, T) — ncq{q^ T) we see that to find the finite temperature scaling we need to change d ^ d— z 
for bosons and d ^ d + z for fermions in universal expressions (j26p and p2p . The same applies to the expression for 
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heat ((25|) . For example, for a slow quench starting at the critical point we find, instead of Eq. ([32]) . 

As before, these scalings are generally valid (to leading order in v) only if the corresponding exponents are less than 
two; otherwise the scaling of riex, Q with v becomes quadratic. 

There is actually an additional potential issue where the scaling (l64t can break. This can happen because the 
integrals over different momentum states q can become divergent at small q. This corresponds to emergence in a new 
non-adiabatic regime of response of the system where the corresponding intensive observables become system-size- 
dependent [l^ . This situation usually happens in low-dimensional bosonic systems where such infrared divergences 
correspond to overpopulation of the low energy modes. While there are specific examples of such behavior for slow 
dynamics in weakly interacting Bose gases 12] and for dynamics in the sine-Gordon model near the massive bosonic 
limit il9|] there is no general understanding of this effect. 



IV. GOING BEYOND CONDENSED MATTER 



Here we focus on some possible extensions of the field of application of slow dynamics to other areas of physics, in 
particular, to cosmology and field (string) theory. The purpose of this section is to show that there are close similarities 
appearing between quantum dynamics near various space-time singularities and quantum critical dynamics. We do not 
intend to give a comprehensive overview of new problems nor present details of derivations. We are not experts in these 
fields. Our purpose is only to show that the applicability of our previous discussion extends well beyond condensed 
matter and atomic physics. Our first example will be from cosmology: we consider an early stage of expansion of 
the universe in d-dimensional space-time and assume that it is described either by a standard de Sitter-type metric 
or by its simplest generalization, which includes a slow-roll parameter. The second example arises mainly in some 
recent studies in string theory, where the question is a propagation of strings on some time-dependent background 
geometries. In particular, these geometries can be singular. This type of singularity is somewhat reminiscent of a 
phase transition, but it is different, and is of a completely geometrical nature. However, it can have an interpretation 
in terms of dynamical systems. We consider the simplest illustrative example - quantum evolution in a so-called 
Milnes universe. In both examples the analogue of a QPT is provided by a singularity in the space-time background 
on which the quantum system is evolved. 



A. Adiabaticity in Cosmology 

A beginning of the Universe can also be considered as a (quantum) phase transition. In modern cosmology one 
usually considers a scalar field (a Higgs field) on some curved background evolving according to the Einstein equation. 
A regime of parameters when the potential energy of that field dominates over the kinetic term leads to the rapid 
blowing up of the Universe, called inflation. Inflation can be modeled by an exponential scale factor in the metric 
and is called the de Sitter epoch. The velocity of inflation is given by a Hubble parameter which is assumed to be a 
constant in the de Sitter model. The natural adiabaticity parameter is then given by the ratio of the Hubble constant 
and the mass of the scalar field. The initial point of the expansion corresponds to the singularity in the metric where 
the scale factor a{t) (see below) is very small. It is therefore natural to look into the cosmological evolution problem 
from the point of view of critical dynamics near the phase transition. 

We consider a spherically-symmetric metric in {d+ l)-dimensional space-time described by a single time-dependent 
scale factor 

ds'^ = -df + a{tydx^, (65) 

where x denotes a d- vector. The Hubble parameter is H = d{t)/a{t), where the "dot" denotes a time derivative. We 
consider a massive scalar field $(x, i}in this geometry, which is minimally coupled to the gravity. The corresponding 
Klein- Gordon equation is given by [33| 

^(x, t) + dH^x, t) - -^\7'^<^{x, t) + ni^^ix, t)=0. (66) 
a (*) 

Introducing the rescaled field (j}{x,t) — a'^/'^$(x, t) and making a Fourier transform we obtain 



(67) 
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where 4'k(t) is a Fourier component of (j){x^t), 
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(68) 



cj^ = + k'^/a^, and we introduced a slow-roll parameter e = H/H^. This parameter is usually considered to be 
small during the inflation stage, thus indicating the smallness of the kinetic energy of a scalar field. 

To proceed with quantization, we expand the Fourier components into a time-dependent creation and annihilation 
operator basis, (pk — ^ ^kit) — i^k{t)dk + '0fc(^)Q'I- The quantization is consistent with the following Klein-Gordon 
scalar product: (01,02) — "i J d'^x{2TT)~'^[(j)l{x,t)4>2{x,t) — (j)i{x,t)(j)2{x,t)]. Because of the time-dependence we have 
to distinguish |in) and |out) vacuum states and corresponding operators a"'°"'. The linear transformation between 
these two bases is given by the Bogoliubov coefficients ak,f3k through ^^"'(i) = ak{t)4'™{t) -\- /3^0)f*(t), whereas the 
relation between vacuum states is given by the squeezed states, |Ofe,in) = |afe|~^ exp(— (/?fe/a^)a™*^a^"^*^)|0/c, out). 
The alternative definition of the Bogoliubov coefficients can be then given via the Klein-Gordon scalar product: 
o-k = (0fc", 0fe"*) and /3fc = — (0™*, 0fe"')- The number of excitations is now naturally defined as 

(Oi„|«r^ar |0i„) = (69) 

The flat de Sitter evolution is defined as a condition that H he a constant, which implies for the scale factor a{t) — 
H^^ exp{Ht). Apparently, in this case e = 0. The solution for the Klein-Gordon equation is then given by (t>^ki^) ~ 
y^n/4:He~'"'^^'Hil\k/[a{t)H]) where n[l\z) is a Hankel function of the first kind. Here 1/ = y/{m/Hy- (d/2)2. 
Using the asymptotics of the Hankel function for large \z\ we can check that this solution indeed describes an |in) 
oscillating state in conformal time r] ~ — exp(— iJt), according to the picture of Ref. [11]. The out solution has the 
following form: 



where Jiu{z) is a Bessel function. Now, the small-argument asymptotics corresponds to the |out) state. Using the 
relation between the Hankel and Bessel functions, 

e-H!:)(.)^ -^^^-(-)-f:^-(-) (71) 
smh(7ri^) 



we obtain the mode-independent Bogoliubov coefficients 

1 



/3 = aexp(-7rz/). (72) 



•^1 — exp(— 27ri/) 
We consider now a regime where [31 

V = H/m < 1 (73) 

which we identify as an adiabatic regime for reasons which will be clear immediately. In the regime of validity of (j73p , 
1/ sa 1/u and one obtains a density distribution of a pair production process at late times (77 — >■ 0) in the form of a 
thermal spectrum 

nex = l/[exp(27r/^^) - 1]. (74) 

The probability of transition between asymptotic in and out vacua is given by the overlap, Pin-s-out — Y\k and 
therefore in our case, the probability per each mode k is given by 

P^^l>out = 1 - exp(-2^/rO, (75) 

which is nothing other than the LZMS transition probability [s^. As we approach the Big Bang singularity the 
parameter v starts to diverge and we expect that the adiabaticity conditions become violated. We thus expect that 
the excitations are created at a much higher rate near the singularity, quite similar to what happens in adiabatic 
dynamics near QCPs. Another possible source of adiabatic/non- adiabatic effects can be found by extending the simple 
de Sitter solution. In the initial stage of inflation, when the potential energy dominates the kinetic term, their ratio 
deflnes a so-called slow-roll parameter which is considered to be small; this parameter is strictly zero in the de Sitter 
case, while it is flnite in more general metrics. This slow-roll parameter can play the role of w, which defines the 
degree of non-adiabaticity in the system. 
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B. Time Evolution in a Singular Space-time 



The idea we put forward in this section can be summarized as foUows: suppose we have a quantum system which 
evolves on a curved space with a time- dependent metric. Suppose that this metric has singularities in some finite 
number of points. It is clear intuitively that the presence of these singularities must inevitably appear in the dynamics 
of the system. In some cases, as we will demonstrate below, nontrivial geometry can mimic nontrivial dynamics, 
similar to dynamics across a QCP. A time-dependent metric of the background geometry plays the role of changing 
the external parameter and thus can induce non-adiabatic effects in quantum dynamics. 

Our simple model here is inspired by recent interest in string theory literature on dynamics of quantum fields, mainly 
of string origin, in time-dependent geometries which contain some singularities. Examples include time-dependent 
orbifolds, null-branes, pp-wave geometries. Big Crunch and Big Rip singularities (see Refs. [40,-^]). The simplest 
illustrative example is the so-called Milne geometry described by the following metric in 1 -I- 1 dimensions: 

ds^ ^ -^dt^ +vH^dx^. (76) 

We consider a quantum scalar field propagating in this metric. Its action is given by 



S =^ I dx I dt\tv\ 



(77) 



where we assume that the evolution starts at i = — cxd and ends at t — +oo. We introduced the scale factor v to 
emphasize its role as the adiabatic parameter. In principle it can be set to unity by appropriately rescaling space-time 
units. The quantum Hamiltonian corresponding to this action is therefore 

^ = 2R / + ^^-'^^'^ + ^ / ^^^^ 

where 11 is the canonically conjugate momentum. 

Apparently, this type of a system can be given a simple condensed-matter interpretation: a massive scalar field 
describes a large variety of one-dimensional phenomena. For example, a Luttinger liquid described by the interaction 
parameter K = \vt\ perturbed by some relevant perturbation (which in the strong-coupling regime of corresponding 
RG can be approximated by a quadratic massive term) could provide a realization of one of the physical models. In 
our case, both the interaction parameter K and the strength of the perturbation depend explicitly on time. From the 
flat-space point of view it is therefore a non-equilibrium model. 

The equations of motion for the model can be put into the form of time-dependent oscillators describing different 
momentum modes: 

<fk+^k{t)^k = 0, (79) 

where ip = ^/tcj) and il.k{t) — m? + [(fc/w)^ + l/4]/t^. The solution can be given in terms of the Bessel functions of the 
order v = ik/v oi the argument mt. 

We note that the situation is rather similar to the previous subsection, except that in the present case the quantum 
evolution goes across the singularity. The presence of the singularity implies that the wave functions of the system 
for t < and for t > must be properly defined and related to each other at the singularity. Several different ways to 
do so have been suggested in the literature: geometric (by going into covering space), operator-analytic (by properly 
regularizing the singularity), dimensional, etc. Without discussing this issue in detail (physically- relevant quantities 
should not depend on the regularization prescription anyway) let us simply state the Bogoliubov coefficients, which 
can be defined again by the matching of in and out states: 

cob{A — ink I v) cos(A) 

■ \,i 1 I \ ' Pfc ^ ^ . , . , , N (80) 

sinh(7rfe/u) sinh(7rfc/i;) 

where the pure phase factor A plays the role of the reflection coefficient and can be defined from the mode-matching 
condition at the singularity, A = aig[Jik/v{nT,t) / J_ikf^{—mt)]. 

The number of particles produced is given by |/3fcp and can be characterized as a thermal distribution with some 
mass-dependent temperature T(rn). For large m, T{m) approaches a constant independent of k and therefore riox ^ 
e~'^™ whereas for small m it is independent of m and n^^ ^ ^ Therefore in the latter case the total number 

of particles N — '^j^n^^ik) scales as \v\. Correspondingly, the in-out transition probability can be put into the 
LZMS form with fc-dependent effective velocity. For large m, P ^ 1 — e^^™ independent of fc, whereas for small m, 
P ~ 1 — g-^^f^/^ independent of m. 
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We therefore conclude that evolution on a singular time-dependent manifold is to some extent equivalent to quantum 
non-equilibrium dynamical systems evolving across a phase transition. In some cases non-equilibrium dynamics of 
quantum system with time-dependent parameters can be modeled by the evolution of a quantum system with time- 
independent parameters on a time-dependent background geometry. This geometry can have singularities which then 
correspond to QCPs in the quantum system. 

V. SUMMARY AND OUTLOOK 

In this chapter we gave an overview of different connections between certain universal equilibrium and non- 
equilibrium properties of continuous quantum phase transitions. It is well known that such transitions are typically 
characterized by singularities in thermodynamic quantities, in particular, in various susceptibilities. We discussed here 
that these singularities lead to the universal nonlinear dynamical response of these observables to various dynamical 
processes near a QCP. 

In particular, we analyzed in detail two possibilities: a sudden spatially uniform quench starting at a QCP; and slow 
passage through a QCP. In the latter case one can also start (end) right at the QCP. In low dimensions such quantities 
as heat (excess energy) generated in the system or density of generated quasi-particles become universal non-analytic 
functions of the adiabatic parameter, which is quench amplitude for sudden quenches, quench rate for slow linear 
quenches, etc. In particular, if the quench \{t) ^ vV then the quasi-particle density scales as riex ~ A 
similar scaling is valid for heat and for other thermodynamic quantities. These scalings are expected to be generically 
valid as long as the corresponding exponents of v remain smaller than two. Otherwise, the leading asymptotic behavior 
becomes quadratic, consistent with standard perturbation theory: ricx ~ v'^ ■ Crossover from non-perturbative to 
perturbative scaling happens precisely at the point where the scaling dimension of the corresponding generalized 
adiabatic susceptibility X2r-(-2(A) (see Eq. (ITUl) ') becomes zero, i.e., when it becomes finite at the QCP. For sudden 
quenches the relevant susceptibility determining nex^ X2('^) is nothing but the fidelity susceptibility; the latter became 
recently an interesting new measure of quantum criticality, independent of the choice of the observable. For slow 
quenches characterized by the exponent r the relevant susceptibilities determining the scaling of the quasiparticle 
density and the probability of exciting the system X2r-i-2(A) is thus a direct generalization of the fidelity susceptibility. 
In particular, the scaling for linear quenches "ex ~ which can be also explained by the Kibble-Zurek 

arguments, [a, [13, [IE], is associated with the scaling dimension of the susceptibility X4(A)- It would be very interesting 
to analyze situations in the future in which many-body excitations do not directly correspond to a fixed number of 
quasi-particles. In those situations one still expects that the scaling for the density of defects will deviate from 
Eq. (|32l) . A possible candidate for this scenario is the non-adiabatic regime, where the density of created excitations 
diverges in the thermodynamic limit [l^ . 

We illustrated how scaling results emerge from adiabatic perturbation theory. This theory has a purely geometric 
interpretation since time can be dropped completely from the analysis (it only implicitly enters through the rate of 
change of the coupling). We discussed a somewhat simplistic situation where only one coupling changes in time. One 
can imagine more general scenarios where several coupling constants change simultaneously. Then the expression for 
the transition amplitude within adiabatic perturbation theory represents a contour integral in the parameter space. In 
this case the expectation values of dynamical quantities will be related to the quantum geometric tensors in the Hilbert 
space of a system. These objects define a structure of the (complex) Riemannian metrics on the space of parameters 
of a system. Evaluated in the ground state, the real part of this metrics is related to the (generalized) fidelity 
susceptibilities, whereas its imaginary part is related to the adiabatic Berry curvature. In adiabatic perturbation 
theory the dynamical phase should then be modified by inclusion of the geometric Berry phase. The Berry curvature 
can diverge close to QCPs. It is then clear that universal features of evolution close to the phase transition advocated 
here should be corrected by specifying the path of the quantum evolution in the parameter space. In particular, 
one can expect additional corrections to the scaling laws coming from the singularities of the Berry phase [sot . We 
expect an interesting interplay between dynamical and geometrical effects in the scaling dependence of quantities in 
the linear-quench regime. Such a possibility to reduce critical dynamics to statics quantities, like quantum geometric 
tensors, looks very intriguing and perhaps requires a closer look beyond adiabatic perturbation theory. 

We also showed that the universal dynamical response can be strongly affected by initial thermal fluctuations and 
that quasi-particle statistics changes the scaling laws. We discussed somewhat simplistic situations where quasi- 
particles are either non-interacting bosons or fermions. In general, critical dynamics at finite temperatures in the 
quantum critical region remains an open problem and it is clear that extra input from equilibrium properties to the 
dynamical response is needed, as compared to the zero temperature case. 

In the last section we gave a brief outlook of connections of critical dynamics with other areas of physics like 
cosmology and string theory. Dynamics near QCPs is qualitatively similar to dynamics near various space-time 
singularities. In some cases non-equilibrium dynamics of a quantum system with time-dependent parameters can be 
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modeled by the evolution of a quantum system with time-independent parameters. However, the background geometry 
of space-time is explicitly time-dependent. Quantum dynamics close to singularities of this dynamical geometry may 
correspond to crossing the phase transition in the static geometry. 

There are many other open questions remaining. The main purpose of this chapter was to shed light on some 
non-equilibrium universal aspects of dynamics near critical points beyond perturbation theory and show their close 
connections to static equilibrium properties. We hope that this chapter will partly stimulate further research in this 
exciting new area. 
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[44] Note that we use dimensionless notations. In physical units the dimensionless combination is m<? j (hH), and thus estabhshes 
a ratio between the Compton wavelength h/mc and the Hubble length c/H. Although H is not a constant, H = H{t), in 
the present cosniological epoch it is iif = 74.2 ± 3.6(krn/s)/Mpc. The adiabaticity condition is thus justified sufficiently far 
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